Constraints on the dark energy using multiple observations : snare of principal 

component analysis 
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We explore snares in determining the equation of state of dark energy (cu) when one uses the 
so-called principal component analysis for multiple observations. We demonstrated drawbacks of 
principal component analysis in an earlier paper. We used the Hubble parameter data generated 
from a fiducial model using the so-called Chevallier-Polarski-Linder parameterization. We extend our 
previous consideration to multiple observations, the Hubble parameter and the luminosity distance. 
We find that the principal component analysis produces the almost constant uj even when a fiducial 
model is a rapidly varying ijj. Thus, resolution of dynamical property of u} through PCA is degraded 
especially when one fits to several observations. 
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The existence of the energy component with negative pressure {i.e. dark energy) besides the matter component 
is one of the alternative elucidations of the current accelerating expansion of the Universe. The lack of enthralling 
theoretical model enforces us to use its equation of state (EOS) cj to classify the different dark energy (DE) models. 
By using the accumulating high precision observational data, the dark energy can be properly investigated. The 
proper paramctrization of uj is prerequisite for fitting the related parameters to data. One of the model-independent 
reconstruction of w is to approximate it by using the piecewise constant bins [1, 2]. It is claimed that one can 
' reconstruct the time dependence of uj and make further model independent studies by using a principal component 
Oh! analysis (PCA) method. This method identifies the directions of data points clustering in the parameter space, 
Q ' and then find a viable dimensional reduction in the parameters with as minimum an information loss as possible 
, [3]. However, it is argued that PCA is sensitive to both the order and the basis set choice [4]. Our conclusion is 
C/3 ' against this argument. We found that PCA produces the almost same results independent of the order of principal 
, components in Ref. [5]. We will also show that PCA gives the almost identical results when we change the basis 

set. One of the main reason for using PCA method is known as to determine whether the dark energy density is 
evolving with time or not. However, it is shown that the time varying ui can be confounded with the incorrect value 
■ of constant one when we use PCA to fit the Hubble parameter H{z) data generated from a fiducial model [5]. We 
adopt the so-called Chevallier-Polarski-Linder (CPL) parameterization oj = ujo + uJa-^ [6, 7] as a fiducial DE model. 
In addition to the Hubble parameter, we use the luminosity distance Dl{z) data generated from a fiducial model to 
investigate the drawbacks of PCA method. 

One divides the redshift range of the survey (z = 0, Zmax) into N bins of not necessarily equal widths Azi {i = 
1, • • • ,N), where J2i — -Zmax- Then, a set of N values of observations of possibly correlated variables can be 
orthogonally transformed into a set of j values of uncorrelated variables so-called principal components without losing 
the original information much [■'>, The dark energy is parameterized in terms of w(z), which is defined to be a 
constant in each redshift bin, with a value uji in ith bin. For the piecewise constant w(z), the energy density of the 
dark energy for z in bin j evolves as 




3(1+Wi) 



where Zi is the lower redshift bound of the ith bin and uji is the fiducial value of the EOS in that bin. Again, we 
emphasize that cj^s have only j — 1 degree of freedom because pde{z) should be equal to pde{z ^ ^) when z = ["]. 
From the above equation (1), this is given by 



Wo 



l+('X:(l+w.)h.[i±^^]) /(ln[l + z,]). (2) 



1=1 

Thus, Uj is determined by other parameters Ui, Zi, and 
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TABLE I: Parameter values of fiducial models are shown in the first two columns. Zi is the uncorrelated bin and there are 
three ujiS at each bin obtained from the minimum fitting by using the generated H{z), Dl{z), and both of them. aiS are 
the l-o" errors of uJiS for corresponding fittings. 
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FIG. 1; Comparison between a fiducial u and ujiS obtained from the different observations, a) The fiducial model is oj = 
— 1.1 + 0.5y^ (dashed) and the obtained values of oJiS from PCA with H data (dotted), Dl data (dot-dashed), and H + Dl 
data (solid). Error bars are obtained from the analysis of -ff + Dl data, b) The fiducial model is oj = — 1.1 
and the obtained values of WiS from PCA with H data (dotted), Dl data (dot-dashed), and H + Dl data (solid) 



0.3 (dashed) 



In vifhat follovifs, we use equally binned both the 40 Hubble parameter H{z) data and the 2000 luminosity distance 
Dl{z) data to cover the redshift range < z < 2 generated from the fiducial models given by CPL parametrization. 
We assume that the errors on both measurements as 5 % and both the present energy density contrast of the matter 
and the present Hubble parameter are exactly known values. We also consider only the flat universe. We 
perform the simple tests to determine the best fit values of tOi and l-tr error which is obtained from the covariance 
matrix. By adding the luminosity distance in the x^ test, we obtain the more confounded results than the one without 
it. This is expected due to the multi-integral of w in it. We show the results for the two different models in table I. In 
the first case, a fiducial model is {ixiq^loo) = (—0.9, 0.5). We perform x^ test with both H{z) and Dl{z) data generated 
from this fiducial model. The fiducial model of the second case is (wq, Wq) = (—1.1, 0.3). We also investigate the PCA 
result of the first model with the different bins. We explain the details of the data with Fig. 1 and Fig. 2 later. We 
choose the uncorrelated bins as Zi = (0.1, 0.4, 0.8, 1.25, 1.6) in the first two cases and Zi = (0.2, 0.5, 1.0, 1.4, 1.6) in the 
third case. As we show before, there are only 4 degree of freedom in this case. Thus, the value of oj^ at z = 1.6 is 
derived from the other w,, z^, and 2:5 values by using Eq. (2). 
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FIG. 2: Comparison of the PCA results for the two different sets of principal components from a fiducial model (ojoji^a) = 
(-1.1,0.5) (dashed), zi = (0.1,0.4,0.8,1.25,1.6) (solid), (0.2,0.5,1.0,1.4,1.6) (dotted). 

We find several snares in PCA when we try to fit both H and Dl data. Firstly, if we exclude the annoying value 
ws, then the dynamical property of w might be detected with PCA from H data as shown in the first fiducial model 
(dashed line) even though the uji values (dotted lines) depart from the values of the fiducial model. However, PCA 
produces the confounded results when we fit uJi to data (dot-dashed lines). The result can be interpreted as 
Lo = —1 for the entire region of z within 1-a error. We reach to the same conclusion when we fit LOi to both H and 
(solid lines). This is shown in the left panel of Fig. 1. The fiducial model is lo = — 1.1 + 0.5 but the result from 
PCA is consistent with a; = — 1. Thus, PCA with multiple observations degrades the resolution of uj. One remark is 
that one might misinterpret that PCA is at least good for fitting H data to reveal the time variation of uj. We showed 
that even the w obtained from the H fitting shows the discrepancy with the fiducial ones [5] . Secondly, PCA produces 
the smaller uji value than the fiducial one when uJa is negative. This is similar to the problem of CPL parameterization 
because it might mislead the property of the fiducial model. It is well known that the negative tOa is harder to be 
detected compared to the positive one [9]. In the right panel of Fig. 1, cuiS increase as z does even though the fiducial 
model UJ = —1.1 — 0.3t^ decreases as z increases. The result of PCA in this case also can be mislead to the constant 

l+z 

u!. Thus, PCA method produces the totally different behavior of uj when uJa is non-negligible independent of the sign 
of it. We can compare this result with one in Ref. [10]. Even though the result in the mentioned reference seems to 
be consistent with the cosmological constant, there still can be the viable time varying DE models which can mimic 
A. This impedes any proper interpretation of the result obtained from PCA. We check that PCA method can give 
the reliable result only when uj is almost constant. 



We show the PCA results from the two different sets of principal components in Fig. 2. We show the results using 
the H + Dl data only. The dashed line is the fiducial model uj ~ —1.1 + O.Sj^ to generate the H and the data. 
The solid lines depict uJiS obtained from binning z as (0.1,0.4,0.8, 1.25, 1.6). The dotted line correspond to the case 
when z is divided by (0.2, 0.5, 1.0, 1.4, 1.6). The results from two different PCs are almost identical. 

Even though; PCA is the most model independent method for probing DE, we demonstrate that PCA method may 
mislead to the property of dark energy independent of the dynamics of the fiducial ut. When one use PCA to several 
observations, the resolution of uj is degraded. The principal component is automatically decided when the data is 
given and our conclusion is independent of PCs. PCA is adequate only when a; is a constant. Thus, we may need to 
check both the model dependent uj parametrization and PCA method to investigate the DE properly. 



4 



Acknowledgments 



This work was supported in part by the National Science Council, Taiwan, ROC under the Grant NSC 98-2112-M- 
001-009-MY3. 



[1] J. Weller and A. Albrecht, Phys. Rev. D 65, 103512 (2002) [arXiv:astro-ph/0106079]. 

[2] D. Huterer and G. Starkman, Phys. Rev. Lett. 90, 031301 (2003) [arXiv:astro-ph/0207517]. 

[3] M. Tegmark, A. N. Taylor, and A. F. Heavens, Astrophys. J. 480, 22 (1997) [arXiv:astro-ph/9603021]. 

[4] T. D. Kitching and A. Amara, Mon. Not. Roy. Astron. Soc.398, 2134 (2009) [arXiv:0905.3383]. 

[5] S. Lee, [arXiv: 1005. 1770]. 

[6] M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10, 213 (2001) [arXiv:gr-qc/0009008]. 
[7] E. V. Linder, Phys. Rev. Lett. 90, 091301 (2003) [arXiv:astro-ph/0208512]. 
[8] K. Pearson, Philosophical. Magazine. 2, 559 (1901). 

[9] L Maor, R. Brustein, J. McMahon, and P. J. Steinhardt, Phys. Rev. D65, 123003 (2002) [arXiv:astro-ph/0112526]. 
[10] P. Serra, A. Cooray, D. E. Holz, A. Melchiorri, S. Pandolfi, and D. Sarkar, Phys. Rev. D 80, 121302 (2009) [arXiv;0908.3186]. 



